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TRANSMISSÃO DO CALOR NUM TUBO EM ESCOAMENTO 
LAMINAR E REGIME NÃO ESTACIONÁRIO 


por JJ. J. DELGADO DOMINGOS 


Prof. do 1.8.T. 
do N. E. E.M. do J. A.C. 


RESUMO 


Apresentam-se as soluções analíticas para o campo das 
temperaturas num fluido em escoamento laminar, para con- 
dição inicial função arbitrária do espaço, temperatura de 
entrada do fluido e temperatura ou fluxo térmico impostos na 


superfície do tubo, funções arbitrárias do espaço e do tempo. 


O método que se apresenta permite utilizar soluções co- 
nhecidas da equação da condução do calor a uma dimensão» 


SYNOPSIS 


This paper presents the analytical solutions for the field 
of temperatures of a fluid in laminar fAlow when the initial 
condition is an arbitrary function of space and the inlet 
temperature of the Auid and the temperature or thermal flux 
imposed on the surface of the tube are arbitrary functions 
of space and time. 

The method here presented allows the utilization of known 


solutions of the one-dimension heat conduction equation, 


0. INTRODUÇÃO 


Não existindo soluções exactas da equação fundamental da transmissão do calor para o escoás 
mento laminar em regime não estacionário, nem valores experimentais que permitam aferir do erro 
em que se incorre baseando a análise no conceito e valores do coeficiente de transmissão por con- 
vecção em regime estacionário, afigura-se de interesse apresentar as soluções intermédias que se 
obtêm admitindo que não existe condução axial e que a velocidade é constante na secção de escoa- 
mento. Estas hipóteses permitem reduzir o problema a um de condução do calor e estender a 
domínio diferente soluções já conhecidas. 

Neste trabalho apresentam-se as soluções formais para o caso geral de a temperatura ou o 
fluxo térmico impostos na superfície do tubo serem funções arbitrárias da coordenada axial e do 
tempo. De igual modo, a velocidade ou a temperatura de entrada do fluido são funções arbitrárizs 
do tempo e do espaço e do tempo, respectivamente. Posteriormente apresentar-se-á uma exemplifica- 
ção numérica das soluções obtidas e far-se-á uma comparação com os valores que resultariam utili- 
zando o método habitual de reduzir a descrição do sistema a uma dimensão e usar valores para o 
coeficiente de transmissão do calor por convecção obtidos em regime estacionário. 
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1. FORMULAÇÃO DO MODELO 
1.0. Admitiremos que: 


a) a secção do escoamento é constante 

b) o sistema possui simetria radial 

c) a velocidade do fluido é positiva e só função do tempo 

d) a condução axial do calor e a dissipação viscosa são desprezáveis 

e) a difusibilidade térmica e a massa específica do fluido são constantes 

f) na superfície cilíndrica que contém o escoamento a temperatura ou o fluxo térmico são fun- 
ções do espaço e do tempo que satisfazem às condições de Dirichelet 

2) a temperatura do fluido à entrada do tubo é uma função do tempo e do raio que satisfaz 
às condições de Dirichelet 


e adoptaremos a seguinte simbologia mais importante : 


T — temperatura 

9) — temperatura do fluido definida pela transformação 3) 

X, x — coordenada axial cuja origem se situa na entrada do tubo e é positiva no sentido do 
escoamento 


r  — coordenada radial 

R -— raio do tubo 

u  — velocidade do fluido 

K  — condutibilidade térmica 
a  — difusibilidade térmica 
t, € — tempo 


Xo — posição dum subsistema em t = 0 

to — instante em que um subsistema entra no tubo 
“ — derivada substancial ou segundo o movimento 
à — Laplaciano. 


1.1 — Atendendo a e) a equação fundamental da transmissão do calor [1] reduz-se a 
1) pc T' = div (K grad T) 
e se a condutibilidade térmica for constante 


2) T = AT 


Se K = K (T), introduzindo a mudança da variável 


T 
3) 9 = — | KdT (*) 
Ks 


ao 


onde K, é a condutibilidade térmica para T = 0, resulta 


4) (5) = K É ie 


(*) Como imediatamente se conclui, se K não é função da temperatura, resulta O = T 
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grad T = e grad €) 


pelo que 2) se reduz a 
5) (5) = A à 


sendo « constante, por hipótese 


1.2 — Atendendo a b) e d) e adoptando uma descrição em coordenadas cilíndricas 5) pode 


escrever-se 
1 
x 


sendo, por c) 


(- , 6) (3 
7) (3) — DO — 0: u é 
Dt dt Jd x 


ou seja, numa descrição de Euler, a equação fundamental da transmissão do calor reduz-se a 


1/90 9 6 O , 1 00 
8) ( ) - + 


u 
ot dx 


1.3 — Por hipótese as condições aos limites são da forma 


9) T= wu (x,t) parar=R t>D0 ex>0 
ou 
dT 
10) do PESO O parar =R t>0 x >0 
e 
11) T'=="Yi tr; É) parar< R t>0 q = 0 


e a condição inicial será 


12) T = O (x, 1) parar<R t=0 x>0 


1.3.1 — As condições anteriores exprimem-se na variável € pela transformação 3) isto é, será: 


9º) B= [PKdT=9 (x, parar=R t>0 ex>0 
Kodo 
o 6) 

10”) — Ko À = q (x, t) parar =R t>0 x >0 
Cr 

; 1 ([" ár E 

11') ara sos KdT=Y(r,t) parar=R t50 dé = O 

Ko Jo 


ap 
12”) D= | 'KdT=P(x,r) parar >R t=0 <> 0 
O u. 


o 
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2, MÉTODO DE SOLUÇÃO (*) 


2.0 Considere-se a fracção de fluido (subsistema) contida entre os planos x, x + dx 
e a superfície cilíndrica r = R. 


du ' 
— == 0, pelo que o subsistema mantém a sua forma geométrica durante 


d 
Por hipótese Gs 
dx 


o movimento no interior do tubo, donde resulta que numa descrição de Lagrange a sua posição fica 
completamente determinada conhecendo a coordenada xo em t=0 (se xo => 0 ou o instante to em 
que o elemento entrou no tubo (se xo <“ 0) ou ainda, dum modo geral, sabendo a posição x que 
ocupa no instante t. 

Mas num sistema de coordenadas cilíndricas fixadas ao subsistema onde fizemos coincidir 
o eixo dos xx com o eixo do escoamento, a evolução das temperaturas no fluido rege-se pela 
equação 6) a qual não é mais que a equação da condução do calor num cilindro cujas bases são iso- 
ladas e onde existe simetria radial. O problema pode pois reduzir-se a um de condução do calor 
desde que no sistema de eixos móveis se introduzam as condições iniciais e aos limites. 


2.1. As condições impostas ao sistema através da superfície do tubo são função da posição 
axial e do tempo. Todavia, em relação ao subsistema, cuja extensão segundo x é infinitesimal, as 
condições impostas em r == R são apenas função do tempo. Com efeito, designemos por f (x, t), sob 
forma genérica, a dependência no espaço e no tempo das condições na superfície do tubo em rela- 
ção a eixos fixos. 


2.1.1. Considerem-se os subsistemas para os quais xo > 0 em t==o, isto é, os subsistemas 
que se encontravam no interior do tubo no instante arbitrário tomado para origem dos tempos. 
No instante t > O estes subsistemas ocupam a região do tubo 


13) <> [ud 


Por outro lado, o subsistema que se encontra em X no instante t encontrava-se em 
14) Xo =X — | udt 


no instante to, e na sua superfície verificava-se portanto nesse instante a condição f (xo, 0). 
Num instante O <. TZ <. t aquele mesmo subsistema ocupa a posição 


t Pr 
-=x— | u dt + | a dt 


at 
15) Es Di | u dt 
A 


ou 


na sua superfície a condição imposta ao instante Z é pois 


16) E (X— [ ua?) 


nd 


* + - . 
(*) O método de solução poderia apresentar-se apenas como uma mudança de variável o que conduziria a uma 
exposição mais condensada, Preferiu-se a via adoptada por ser fisicamente mais intuitiva. 
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2.1.2. Considerem-se agora os subsistemas que entraram no tubo apenas em to > O. Estes 
ocuparão no instante t a região do tubo 


17) «< [ud 


e um subsistema que no instante t se encontra em X entrou no tubo no instante to dado implici- 
tamente por 


t 

18) Xi udt 
to 

Façamos 

19) Cemst-to>oO 


com o que introduzimos uma nova origem dos tempos para cada subsistema cuja posição 
satisfaz a 17). 
Em é =0 o subsistema entra no tubo e a condição na superfície exprime-se por f(O, to). 
Num instante 0O<. T< É encontra-se em 


t+to t 
20) «= [ udt=X— [ udt 
to t=—T 


sujeito à condição à superfície 
t-Hto 
21) ECfo udtit+to) 
to 


2.2. As expressões 16) ou 21) são, relativamente ao subsistema, apenas função de Z e tradu- 
zem, no sistema de eixos móveis, as condições à superfície expressas em eixos fixos pela forma 
genérica f(x, t). 

X e t são parâmetros que identificam o subsistema alternativamente a xo e to, como já se referiu 


2.3. Expressas as condições aos limites a formulação matemática exige o conhecimento das 
condições iniciais para cada subsistema. 


2.3.1. Para os subsistemas cuja posição satisfas à expressão 13) a condição inicial é, por 14). 
t 

22) T=E(X- [ udt, ) = Hi (X t,1) 
o 


a qual, pela transformação 3) se exprime no variável O por 


p 
23) 9==[ iKdt= 0(X,t,1) 
Ko Jo 


2.3.2. Para os subsistemas que em t ocupam o domínio defenido por 17) o instante inicial 
corresponde a Z == 0. Atendendo a 19) e 11) a condição inicial para o subsistema que no instante t 
se encontra em X será pois 


24) T= F (r, to) para t=0 r<R 


ou, na variável €) 
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(et, 
25) A = E) 4 "| KdT para 7Z=0 


sendo to dado por 18). 


2.4. Expressas as condições iniciais e aos limites para qualquer subsistema, vejamos segui- 
damente a solução do problema para dois casos particulares de condições sobre a superfície. 


3. FLUXO TÉRMICO IMPOSTO NA SUPERFÍCIE DO TUBO 


+ . . O +. t 
3.1. Em relação a eixos móveis e para os subsistemas que ocupam o domínio X > f E udt o 


problema formula-se matemáticamente do seguinte modo: 
Determinar O(r,z) tal que 


1 2 (5) 1 dt 
Agra O qua f<R T>O 
& o Tr r dr 
dt) 
— Ko = q (X — udt, Z) 
7, r z 
p=R T>0 
e q (X, t, Z) 
Ea É tm (X, t, 7) Red 
Raso 
= Oo (X, t, 1) IR T=0 
como facilmente se verifica 
27) (3) Tue (O, + (3, 
em que €, é tal que 
1 J2 (3) d (3) ” 
tt (4 — = 1 À r «E R t A O 
x dr? r dr 
3 E) 
28) TD =] r=R 7>0 
dr 
(4 — o (X, E: r) LS R t=0 
e (6 tal, que 
92 6) ) O 
E Pis: Tia À E<CR T>O0 
d Ta r dr 
d (Cy 
29) = Ko ri (X, t, Z) tr=R E%t.5bO 
r 
(o, = 0 r << R t=0 


sendo, pelo teorema de Duhamel |1). 


70 


TE UNICA Ns 352 


“e | 
30) Ds (X, t, F, t) es | O, (r, T — 4) Ed (X, t, ) d À + O; (r, Í) q: (X, t, 0) 
. [8] 


em que (); satisfaz a 


1 o? Es 1 od Os 
2 Of = RE 1 r<R T>O0 
á dr + r dr 
d O 
31) Ko 5 7 r=R 11>O 
r 
O = L< R t=o0 


3.1.1. A solução explícita de 28) depende da forma de O,. Por hipótese Oo satisfaz às con- 
dições de Dirichelet pelo que a solução de 28) é no caso geral (2). 


32) O (X, tr, 7) = al Cn e Jo (An 1) 


sendo 


2 R 
33 Cn o X, t: o À d r 
) R)] ) a ( r) Jo Un r) 


CRE O RIA Qu 
onde ), são as raizes de 
34) H(tR)=o0 


A solução de 31) obtém-se pelo mesmo método e é [3] 


o oc -infat » (à 
35) O (r, 7) = ei a : E ad En º Jo Cat) 
Ko R Ko 2 me 4 n== An? Rº Jo (In R) 


sendo os à, as raizes de 34). 
Substituindo 35) em 30) obter-se-á (Os (X, t, r, Z). Por 27) obteremos & (X, t, r, ZT) que será a 


temperatura em função de r no instante Z para o subsistema identificado por se encontrar em X no 
instante t. 


3.1.2 — As soluções anteriores, onde X e t são parâmetros, dão a temperatura no subsistema 
em função do raio num instante Z. Como porém o subsistema se encontra na coordenada axial X 
no instante t, se fizermos Z ==t na expressão de O (X, t, r, Z) obteremos O (X, t, r) que será a 
temperatura do fluido expressa em relação aos eixos fixos e em função do tempo no domínio defi- 


nido por 13). 
3.2. Para os subsistemas que no instante t ocupam o domínio 


t 
36) X< [ud 


o problema formula-se, matemáticamente, de modo análogo ao anterior: 
Determinar &) (r,z) tal que 


1 2 (5) (5) 
O = o a TS R t> 0 
a dr * r dr 
o E) ttto 
— Ko =a(/, udt,Z+to) r=R t>O 
Or to 
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37) 


ca q tt; to) 
1 “MW (r,to) 
O) = / K dt= O (r to) 
Ro « O 
sendo, novamente 
38) O) — O, - (3; 


em que (O, é dado pela expressão 32) sendo, todavia, os cn dados por 


2 R 
9 1 í-———>——— — So , to o ka d 
a "O R Do? (Un R) + JH? (Mn R)] |, Op No Jo Pan Es 


6h; (Z, r, to) obtem-se de (3 dado por 35) usando o teorema de Duhamel 


T 
40) O; (r, Z, to) = f Os (r, 7 —4) = (A, to) di + Os (1,7) qr (0, to) 


o 


3.2.1 Tal como anteriormente, () exprime a temperatura no subsistema que entrou no tubo 
em t == to em função de r e Z. Relativamente a eixos fixos () exprime a temperatura na secção de 
coordenada 


41) x=[ udt 


to 


no instante t = Z + to, em função de r. 
Como u = u (t) é conhecido, de 41) poderá determinar-se 


42) to = to (X, t) 


que substituído em & (r, Z, to) exprimirá O) em relação aos eixos fixos. 


4. TEMPERATURA IMPOSTA NA SUPERFÍCIE DO TUBO 


et 
4.2 — Para os subsistemas que ocupam o domínio x > ) udto problema formula-se mate- 
(o) 


máticamente como 
Determinar O (r, Z) tal que 


E a J2 E) 1 dO 
x or 1 dO 


at 


O = (X—|] ud) =u(Xt?) r=R => o 
44) ão 


“L 
0= — | K dt = Oo (X,t, r) r< R €t= 0 


Ko é o 
e €) obtém-se, análogamente aos casos anteriores, de 
45) O) = (6 + (O, 


4.1.1. — Oy satisfaz à equacão diferencial 43), à condição do limite de a temperatura ser zero 
em r==R e à condição inicial 
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TECNICA XIV 


46) O = 95 (X, t, 1) 
sendo a solução 
oo -in at 


47) O = 2 oeê Jo (An r) 


Os coeficientes do desenvolvimento em série, cn, determinam-se pela expressão 33) e os valores 2n 
são as raizes da equação 


48) Jo (An R) = 0 
41.2. — €s obtém-se pelo teorema de Duhamel 


ag at 
49) ("es (rt) e 0,4) 40 (1,2) 4X [ ud, 0) 


onde €; (r, Z) é a função que satisfaz à equação diferencial 43), à condição ao limite de O; = 1 em 
r==R e à condição inicial de a temperatura ser zero parar < R. 


(6); é dado por [2): 


- — In* a % 
e 2 Ji (4n R) . e 
so Meat tt 
s (1,2) n=1 Rd [Ji (MR) + Jt (An R)] de Vi) 


4.1.3. Tal como em 3.1.2 a descrição da temperatura do fluido em relação a eixos fixos 
obtém-se fazendo t == na expressão de O (X,t,r, 2). 


t 
4,2. Para os subsistemas que ocupam o domínio x < 7 udt, € obtém-se relativamente a 
o 


eixos móveis pela solução de 


92 6) 1 960 


—-—— 


2 
—— (5) — 
sa) x dr? r dr 


CM, EB 


Tito 
e=2([ udt, 7 + to) = g2(Z,to) p= Ri T>0 
to 


52) 


Wi (r, to ? 
==)" KdT=" (r, to) t<R tT=0 
9 o 


sendo ty O instante em que o subsistema entra no tubo, o qual será dado implicitamente por 


t 
53) x=[ udt 


to 
se o subsistema se encontrar em X no instante t. 
4.2.1. O obtém-se também da soma: 
54) 9= 0 + 0; 
em que O; satisfaz 51) e à condição ao limite 


O, = 0 para r=R e à condição inicial 
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O = O(r,to) parar CR. O (r,to,Z) é dado pela expressão geral 47) em que os 4, são raízes 
de 48) e os cy são determinados da expressão 33) substituindo Co por “F (r,to). 


4.2.2. Os (r,7,to) obtém-se de 


t 
55) = "ney LL amd 
o dz 


sendo (); dado por 50) 


4.2.3. Identicamente ao exposto em 3.2.1 a expressão de €) relativamente a eixos fixos e a t 
obtém-se exprimindo to em função de X pela relação 41) e substituindo em C(r,to, Z), T por 
Cc=t—to. 
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TRANSMISSÕES EM ONDAS DECAMÉTRICAS 


PREVISÃO DA FREQUÊNCIA A UTILIZAR 


RESUMO 


O processo corrente de previsão das frequências a utili- 
sar nos sistemas de comunicação em ondas decamétricas, 
consiste exclusivamente na determinação do valor mediano 
da frequência máxima utilizável e de uma frequência óptima 
de tráfego que para a camada Fa se faz igual a 85º da- 
quela mediana. No presente artigo, expõe-se um método de 
cálculo que faz intervir a distribuição estatística dos parã- 
metros que determinam as condições de propagação, cujos 
resultados, expressos na probabilidade de que a comunica- 
ção se faça com êxito, são de maior confiança, e facilitam 
a escolha correcta da frequência a utilizar. 


1. INTRODUÇÃO 


A atenuação das ondas electromagnéticas que 
se propagam por reflexão nas camadas ionizadas 
da atmosfera (caso das ondas decamétricas), está 
sujeita a variações regulares que, para um dado 
circuito radioeléctrico, são quase exclusivamente 
determinadas por três parâmetros: frequência, 
hora da transmissão e índice de actividade solar. 
Se se considerar a variação das condições de pro- 
pagação ao longo de um mês, período usual nas 
previsões de frequência, podem tomar-se como 
variáveis, apenas a frequência e a hora da trans- 
missão, e exprimir as condições de propagação 
em função destas duas grandezas, procedimento 
que é universalmente seguido. 

A estas variações regulares sobrepõem-se ou- 
tras de ocorrência aleatória, que resultam das 
alterações irregulares das características da ionos- 
fera (densidade de ionização, altura, perfil e forma 
das camadas) e ainda da existência simultânea 
de diversos modos de propagação. O ruído, de 
que depende a qualidade da recepção, está su- 
jeito às variações anteriores e ainda às alterações 
casuais do grande número de fontes que o geram. 
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SINOPSYS 

The process currently used in the forecast of the frequen- 
cies to be employed in broadeasting systems using decame- 
tric waves, consists exclusively in determining the median 
value of the Maximum Usable Frequency ( MUF) and a 


frequency called optimum traffic frequency which, for the 


layer Fo, is made equal to 850/, of that median value, In 
this article a method of calculation that deals with the sta- 
tistical distribution of the parameters that determine pre- 
pagation conditions is exposed, The results of this method 
expressed in the probability of sucess of the communication, 
are more secure than usual in classical methods and they 
facilitate the choice of the correct frequency to be used, 


Estas variações irregulares do sinal e do ruído 
têm influência na qualidade do serviço de comu- 
nicação, mas dada a sua natureza não podem ser 
previstas. 

Trabalhos relativamente recentes do C.C.IR. e 
de outros Organismos que se dedicam ao es- 
tudo da ionosfera, mostraram que as variações 
irregulares do sinal e do ruído podem ser des- 
critas com boa aproximação por leis estatísticas 
simples, o que permite tomá-las em considera- 
ção na previsão das condições de recepção das 
emissões em ondas decamétricas. 

O método de previsão que vamos expor con- 
siste em calcular a probabilidade de que a co- 
municação satisfaça às condições desejadas para 
cada par de valores: frequência — hora da emis- 
são. Efectuando o cálculo para uma rede de pon- 
tos conveniente, podem traçar-se as curvas de 
probabilidade constante num sistema de coordena- 
das horas — frequência ; obtém-se assim um dia- 
grama em que são considerados todos os ele- 
mentos que condicionam a qualidade da recep- 
ção, e cuja interpretação é imediata. 

Um diagrama deste tipo é apresentado por 
Greenberg na publicação «Optimum H.F. Predic- 
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tion», [2], ma elaboração do qual são tomadas 
em consideração as distribuições de dia para dia 
do sinal, do ruído, e da frequência máxima utili- 
zável. 


2. PROBABILIDADE DE SERVIÇO 


Nas transmissões em espaço livre a potência 
do sinal captada pela antena receptora, s, é dada 
pela expressão 

3. 43 


= (5) t Br Br (1) 


" 


em que 


À — comprimento de onda 

d — distância entre as antenas emissora e re- 
ceptora 

t — potência radiada 

Br € Br 7 
isotrópica das antenas emissora e receptora na di- 
recção em que se faz a transmissão. 

Em unidades logarítmicas (*) 


) 


ganho relativamente à antena 


Nas comuniçações em ondas decamétricas en- 
tre dois pontos da superfície terrestre intervêm a 
atenuação ionostérica A, e a atenuação por reflec- 
ção no terreno Ar. Assim teremos 


47% 
95=T-20log, are + Gr + Gr -Aj— Ar (2) 


Nesta expressão a atenuação ionosférica A, é 
uma grandeza com carácter estatístico como 
apontámos em 1. As respectivas variações de dia 
para dia ao longo de um mês seguem com boa 
aproximação a distribuição de Gauss com um 
desvio padrão de 5 dB. 

Sendo constantes os restantes termos da ex- 
pressão anterior, S seguirá também a lei normal 
com o mesmo desvio padrão, e mediana 


E 4 x 
3 = Te 20 loga + G.+ GrR— A, — Ar (3) 


em que A, se calcula por uma expressão semi- 
-empírica devida a Haydon e Laitinen [3]. 


- nn mememe m— 


(*) Representaremos por maiúsculas as grandezas ex- 
pressas em dB, 
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Designando a potência de ruído captada pela 
antena receptora por N, a relação sinal-ruído R 


terá o valor 
R=S—N (4) 


em que N, tal como S, é uma grandeza com dis- 
tribuição normal, cuja variância é função da fre- 
quência, da hora da transmissão e da estação do 
ano. Considerando as variáveis S e N incorrela- 
cionadas, a diferença S — N segue uma lei nor- 
mal com um valor mediano 


R=S—N (5) 
e uma variância 
or? = 05? + on? (6) 


Sendo Ro a relação sinal-ruído mínima, que a 
qualidade de recepção pretendida impõe, a dife- 
rença 


R — Ro (7) 


determina a fracção do tempo total em que a 
relação sinal ruído Ro é igualada ou ultrapas- 
sada. Esta fracção é designada duração de ser- 
viço, nas publicações do C.C.I.R. Será represen- 
tada por à (fig. 1). 


o 
ta 
> 
E4 
(a) 
“ 
(o 
La 
o 
s 
Õ | 
ns 
es e amet Saes 
: =. 
ie 
o Ro [o R 
Fig. 1 


Como a distribuição de R é normal, à satisfaz 
a relação 


É 1 [ + 00 
== FE Cx aa 
V2m GR P 


d (e) 


Introduzindo nesta expressão um parâmetro m 


e uma variável x definidos por 
R—R 


R—-R6=mo (8) x = (8") 
SR 
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resulta 


am 3 
PEDE. | exp (— =) ds (9) 
VIE da és 2 


Designando por m, o valor de m que corres- 
ponde à duração de serviço desejada, e por S, o 
valor mediano da potência do sinal que garante 
esta duração de serviço, de (5) (7) e (8) obtém-se 


So = N + Ro + mo 9R (10) 


Se se pretender operar com a potência do si- 
nal relativa à portadora, como é usual, a relação 
Ro que intervém nesta expressão tem o signifi- 
cado de relação (em dB) entre a potência da 
portadora e a potência do ruído captadas pela 
antena receptora. Ro depende da relação sinal- 
-ruído desejada em audio-frequência, que carac- 
teriza a qualidade da recepção; esta relação pode 
fixar-se a partir dos resultados de ensaios de 
perceptibilidade (palavra), e apreciação subjec- 
tiva (música), efectuados com grande número 
de indivíduos. A relação Ro é função da relação 
sinal-ruído em audiofrequência e do factor de 
modulação (V. apêndice). 

Em geral, o valor da potência necessária à sa- 
tisfação da qualidade de serviço desejada, dado 
pela expressão (10) será diferente do valor me- 
diano da potência captada pela antena receptora, 
dado por (3). A diferença entre estes dois valo- 
res 

D=S-S=S-N-R -moor (11) 
dá-nos uma indicação do êxito da transmissão. 
Contudo, quer S quer So estão sujeitas a desvios 
casuais que lhes conferem o carácter de variá- 
veis estatísticas. Assim 


— S está afectada das incertezas dos elementos 
que intervêm no cálculo da expressão (3). 
No caso da radiodifusão refere-se em parti- 
cular a incerteza relativa a Gr. 


— N apresenta desvios de ano para ano e de 
local para local. 

— Ro depende do ouvinte, e das características 
do programa (caso da radiodifusão). 

— qr também está afectada de incerteza. As 
leis de distribuição das variações de dia 
para dia do sinal e do ruído são leis mê- 
dias, relativamente às quais se observam 
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desvios de local para local e de ano para 
ano. 


Para tornar mais clara a escrita da equação (11) 
substituiremos S por S -+- 4Se Npor N + AN 
passando S e N a terem o significado de média 
das médias. Deste modo, 


D=S-+AS-(N+AHAN)-Ro-—- mosr (12) 


Com aproximação suficiente pode admitir-se 
que as grandezas que intervêm nesta expressão 
seguem a lei de Gauss o que significa que, com 
melhor aproximação, D terá a mesma distribui- 
ção. 

Tratando-se de variáveis sem correlação ob- 
tém-se 


D =S-N-Ro—mo?R (13) 


Tp E Tg” + TN + Ro. + mo” TOR (14) 


em que 252, 05º, ON, Ros Cop: SãO as variân- 


cias respectivamente de D, AS, AN, Ro € GR. 


Desde que seja possível determinar os valores 
de todas as grandezas que intervêm nas duas 
expressões anteriores, a probabilidade de que a 
relação sinal-ruído necessária à transmissão seja 
ultrapassada durante a fracção do período de 
previsão, calcula-se imediatamente a partir da 
distribuição normal: é a probabilidade de que a 
diferença D seja positiva. Representando-a por 
Ps, tem-se então 


1 Rs” (15) 
RE I exp 
] 


V2% 9%. 


(D-DP 


“20? 


Ps = jap 


cujo valor se obterá numa tabela da lel de Gal- 
> E) 


ton passando para a variável reduzida 


A probabilidade Ps é designada probabilidade 
de serviço no Rapport 322 do C.CIR. 


+ + * 


Passamos a seguir ao cálculo dos termos do 2.º 
membro das expressões (13) e (14). 

S, Ro e mo já anteriormente foram referidas, 
e 7x é dado no Rapport 322 do C.CI.R. [1] 
pelo que somente trataremos das restantes gran- 
dezas : 
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N — Em tudo o que atrás dissemos, ignorámos 
o ruído interno do receptor, porquanto supuse- 
mos idênticas a relação sinal-ruído nos terminais 
da antena receptora e a mesma relação imediata- 
mente antes do detector. 

O ruído introduzido pelo receptor pode ser 
tomado em consideração, tratando o conjunto 
antena — linha de transmissão — receptor, como 
uma cadeia de quadripolos alimentada por ruído 
térmico e com um factor de ruído equivalente 


fea — mea * E Fa 


im 


em que To é a temperatura dos quadripolos, Ta 
a temperatura equivalente da antena, e fy o fac- 
tor de ruído da cadeia quando a temperatura da 
fonte é igual à temperatura dos quadripolos. 
Tudo se passa como se não existisse ruído in- 
terno, desde que se tome para potência captada 
de ruído 


N = 10 log, (K To B feq) 


em que K é a constante de Boltzman e B a lar- 
gura de banda do amplificador de frequência in- 
termédia do receptor. 


- F = ER 
Na faixa das ondas decamétricas — » fy — 1, 
E) 
o que permite desprezar o ruido interno. 


No Rapport 322 do C.C.I.R. são dados os va- 
lores de Ta sendo To=-278'K. 


a) 
3R e Cor — À média e a variância de cr são fun- 
ções da média e da variância de cs e 7N através 
da relação (6). Dada a sua natureza, nenhuma 
destas grandezas pode satisfazer a lei de distri- 
buição normal, mas a dispersão dos valores ex- 
perimentais é suficientemente pequena para per- 
mitir substituir a lei real pela lei normal no in- 
tervalo entre os decis. A concentração de valo- 
res torna ainda possível o desprezo dos termos 
de grau superior ao primeiro no desenvolvi- 
mento de 7R em série de Taylor em torno do 
ponto (75, 9N). 
Deste modo obtêm-se de (6) 


Le] a 


E E E R Es 
aR=Vas + ON + pasa == (35 — 05 + 
É e Nr Ru 
Vas +an 
EN o 
+ ==". (uh = 8) 
| a hr) 
Vos + an 
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de que resulta 


or CoR=05 Cos + GN Toy 
e = E) ; a A “ d g 5 aaa 
em que 75; e 73; — variâncias de 05 e GN 
podem obter-se a partir das publicações do €.€. 
LR, 


os e cp,-- Os valores destas grandezas terão 
que ser estimados podendo adoptar-se 1 dB para 
a primeira e 2 dB para a segunda. 


3. PROBABILIDADE DE ÊXITO 


A probabilidade de que uma comunicação via 
ionosfera se faça com êxito depende, por um 
lado da possibilidade de que seja atingida a rela- 
ção sinal-ruído necessária durante a fracção do 
periodo de previsão imposta pela qualidade de 
serviço desejada (probabilidade de serviço definida 
em 2.), e por outro, da probabilidade de que à 
frequência de tráfego a ionosfera apresente ca- 
racterísticas refiectoras. Designando por fo a fre- 
quência de comunicação e por É a frequência 
máxima utilizável, a probabilidade de êxito será 
a probabilidade de que 


D>0 e £f-fo >0 


Segundo o National Bureau of Standards, as 
variações de dia para dia durante um mês, da 
frequência máxima utilizável da região F» seguem 
aproximadamente a lei normal com um desvio 
padrão c;=:0,117 É sendo f a média mensal. 

Visto que D e f são variáveis sem correlação 
a probabilidade de êxito tem o valor 


P= Ps x P (18) 


em que Ps é dada por (15), e 


1 | OO 
= exp 


o (E— f)? 
20F 


fo 


de a refiexão se verificar na camada E, o que 
acontece por vezes em ondas decamétricas, será 
P;=1 para É > fo. 

A relação (18) permite determinar, para cada 
par de valores frequência — hora da transmissão, 
a correspondente probabilidade de êxito. Traçando 
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as curvas de probabilidade constante num sis- 
tema tendo por coordenadas aquelas variáveis, 
obtém-se um diagrama como o da fig. 2, que 
torna simples a escolha correcta da frequência a 
utilizar. 

Este método de apresentação das previsões da 
frequência a utilizar tem sobre o clássico (MUEF, 
FOT, LUF) as seguintes vantagens: 

— Os resultados são apresentados numa forma 
mais explícita e condensada. Num único diagrama 
encontram-se todos os elementos que interessam 
à escolha da frequência de tráfego em função da 
hora da emissão, expressos num parâmetro de 
interpretação imediata ; 

30 


22 

20 

“€ I8 
r 
3 

o 16 
“UV 
c 
«o 

T I4 
Prel 

12 

To) 

8 

6 

4 
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Ber es ng 
DRE ARRE 


— é mais informativo porquanto no método 
clássico apenas se determinam limites fora dos 
quais se não deve comunicar; 

— a previsão é de maior confiança porque se 
entra em consideração com as incertezas de to- 
das as variáveis intervenientes e a respectiva dis- 
tribuição estatística, 


Tem contudo o inconveniente de os cálculos 
serem bastante mais extensos, justificando-se a 
suá adopção na hipótese de se dispor de meios 
de cálculo automático, para o que exige um pro- 
grama mais simples do que a determinação da 
frequência útil mínima. 


T. MG. 


Fig. 2 — Probabilidade de êxito de uma comunicação Lisboa - Lourenço Marques em Março de 1965, 
utilizando um emissor de 100 kW e uma antena com um ganho de 20 dB. Adoptou-se para R, 
o valor de 26 dB. 
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4. APLICAÇÃO 


O diagrama da fig. 2 é o resultado da aplica- 
ção ao circuito radioeléctrico Lisboa-Lourenço 
Marques do método de previsão exposto em 2. e 3. 

A observação do diagrama mostra-nos a varia- 
ção da frequência óptima de tráfego acompa- 
nhando a variação da densidade de ionização da 
região reflectora, pondo ainda em evidência a 
impossibilidade de comunicação com êxito entre 
as 03 e 07, conclusão que pode não ser aparente 
num diagrama MUF — FOT — LUF. 

Na fig. 3 ilustramos um outro modo de apre- 
sentação das previsões, útil para a radiodifusão 
à qual estão reservadas determinadas faixas de 
largura estreita, dentro das quais a variação da 
probabilidade de êxito é insignificante. 


p 


LO —-— 
| Mm. 


| 


| 


20 o 4 8 
HORAS T.M.G. 
Fig. 3 
APÊNDICE 
Sejam 
e(t)= SI1-+ucosmst] cos mot 
A (1) 


n (t) = na (t) coswot + np (t) senwot 
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respectivamente as tensões do sinal e do ruído 
à saída do amplificador FI do receptor; wo é a 
frequência central deste amplificador, » é o fac- 
tor de modulação e na (t) e mp (t) duas funções 
sem correlação de variação lenta e com o mesmo 
valor eficaz (Nes. Das relações anteriores con- 
clui-se que a envolvente da tensão de radiofre- 
quência no detector é dada por 


u(t) = V[% + u &ocos os (t) + na (D]! + np 2 (t) 
(2) 


de que se obtém 


E PET R 1 ve 9 
a (3) 


Desenvolvendo (2) em série e limitando o de- 


senvolvimento aos termos de 2.º grau ( - >10) 
Net 
9 
ST E? 1 Ne 
Lo =d ES (4) 
2 “So 


por ser ma(t) = np(t) = 0. De (3) e (4) con- 
clui-se que a relação sinal-ruído em audiofre- 
quência, considerando detecção de pico, tem o 
valor 
28.2 
Rm “o 
p= ——— 
E ao (5) 
Designando por r a relação entre as potências 
da portadora e do ruído antes da detecção de 
(1) e (5) obtém-se 


2 


"o 


15 
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MECÂNICA DOS SOLOS EM REGIME DINÂMICO 
VIBRAÇÃO DE UM BLOCO ASSENTE SOBRE UM MEIO ELÁSTICO 


(Elementos) 


RESUMO 


Caracterizam-se as propriedades elásticas dos solos, o 
dominio de aplicação das leis de Hooke e definem-se, anali- 
tica e experimentalmente, para diversos tipos de solo, os pa- 
râmetros que caracterizam as forças de restituição dos solos, 
sujeitos a vibrações, 

Apresentam-se em seguida aspectos da Teoria das Vi- 
brações, particularizada para o problema do bloco rêgido 
assente sobre um meio elástico e sujeito q vibrações, em re- 
gime livre ou forçado, de diversos tipos, 


O presente artigo constitui o primeiro de uma série 
na qual serão abordados os principais problemas en- 
volvidos no projecto de fundações de máquinas alterna- 
tivas, de percussão e rotativas. 

A resolução daqueles problemas exige conheci- 
mentos em 4 domínios: 

— Mecânica dos Solos, especialmente no que se 
refere a características elásticas dos solos e propagação 
de vibrações; 

— Teoria das Vibrações; 

— Elementos de Máquinas, no que se refere à de- 
terminação das solicitações devidas ao seu funciona- 
mento e transmitidas às fundações; 

— Teoria das Estruturas. 

Da diversidade dos domínios indicados ressalta cla- 
ramente a razão pela qual estes problemas não são 
normalmente encarados na sua verdadeira dimensão e 
se tenta torneá-los reduzindo-os a cálculos estáticos 
correntes, através de coeficientes dinâmicos de valor 
mais ou menos elevado, o que de forma nenhuma asse- 
gura um comportamento eficaz das fundações. 

Daí resultam graves inconvenientes para o funcio- 
namento das máquinas, com sensível redução na sua 
vida útil e acréscimo desnecessário no custo das fun- 
dações cujo sobredimensionamento arbitrário não im- 
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SYNOPSIS 


Elastie properties of soils and the field of application 
of Hooke's laws are here characterized and the restitution 
forces of soils subject to vibrations are analytically and 
experimentaly defined for several types of soils, 

Some aspects of the Theory of Vibrations particula- 
rized to the problem of a rigid block resting upon an elastic 
medium and subject ed to various kinds of transient or 
steady vibrations are then presented. 


pede muitas vezes a propagação de vibrações de am- 
plitude exagerada com a consequente danificação em 
estruturas próximas ou até afastadas muitas dezenas de 
metros. 

Ainda que a análise detalhada dos problemas apon- 
tados exceda largamente o âmbito da presente série 
de artigos, julga-se que esta poderá contribuir para o 
seu esclarecimento, ajudando a formar uma «sensibili- 
dade técnica» em relação ao comportamento das es- 
truturas actuadas por solicitações em regime dinã- 
mico. 


1. CARACTERÍSTICAS ELÁSTICAS 
DOS SOLOS 


1.1 — Aplicação da lei de Hooke 


A lei de Hook aplica-se a sólidos homogéneos 
e isótropos. Relaciona, num ponto qualquer do 
sólido, o estado de tensão com o estado de de- 
formação. Assim, pode exprimir-se pelo seguinte 
conjunto de equações: 
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ox =h6+2U ty 


iy = 4e+ 24 E y 

G =he+24p 

e = +Hiy+tos (1) 
Cxy = 2º Yxy 

Cyz = 2W Yxz 

Cyz = 2H Jxz 


designando 7x, 7y, 7 as componentes normais 
de tensão, Txy, Txz , Tzy as componentes tangen- 
ciais de tensão, “x, Ey, :z as componentes nor- 
mais unitárias da deformação, xy, Yxz, Yyz OS 
deslocamentos angulares unitários, ou desliza- 
mentos, e 4 e as constantes de Lamé, relacio- 
nadas com os módulos de elasticidade longitu- 
dinal E, transversal, G e com o coeficiente de 
Poisson, v, pelas equações: 


ses 

1-+-v (1 —2u) 
(2) 

ap 

2(1 + 2) 


No caso de sólidos anisótropos, são neces- 
sárias mais que 2 constantes independentes para 
relacionar os estados de tensão e deformação. 

A lei de Hooke é válida para o caso de não 
existirem tensões internas iniciais, ou de estas 
serem muito pequenas em relação às tensões apli- 
cadas. Neste caso tanto o módulo de elasticidade 
longitudinal como o coeficiente de Poisson são 
constantes. 

Nos solos correntes existem sempre tensões 
iniciais, devidas às forças de coesão, atrito interno 
e capilaridade. Em certos casos estas tensões ini- 
ciais são de valor muito elevado, e a lei de Hook 
é totalmente ineficaz para traduzir as relações 
entre o estado de tensão e de deformação. É 
este, nomeadamente o caso do interior da Terra, 
onde se verificam tensões consideráveis devidas 
à atracção de grandes massas. 

Mesmo nos casos em que as tensões iniciais 
são pequenas, elas existem sempre. 

Do mesmo modo, o módulo de Young não é 
uma constante, nos solos, mas depende do es- 
tado de tensão. 

As figs. 1 e 2 explicitam respectivamente a 
existência de tensões iniciais e a variação do mó- 
dulo de Young: 


82 


% 


dy 
Fig. 1 — Relação entre a tensão tangen- 
cial e o deslizamento, em solos correntes 


0z 


Ez 


Fig. 2 — Relação entre a tensão normal 
e a extensão, em solos correntes 


Deduz-se facilmente, a partir dos dois fenó- 
menos indicados, que o valor do coeficiente de 
Poisson não é constante, mas varia também com 
a tensão normal. 

Por outro lado, a deformação de um solo de- 
pende do tempo de aplicação do estado de ten- 
são, aumentando com aquele. 

Desde que se tomem em conta os factos apon- 
tados, é-se levado a uma teoria de elasticidade 
não linear, expressa por equações diferenciais 
não lineares, cujas soluções, mesmo nos proble- 
mas mais simples, são extremamente difíceis de 
obter. É esta a razão porque toda a análise dos 
solos é baseada nas leis de Hooke, ainda que 
isso constitua uma aproximação considerável. 


1.2 — Valores do módulo de Young e coeficiente 
de Poisson 


O coeficiente de Poisson é sempre menor nas 
areias que nas argilas. Ainda que os resultados 
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experimentais sejam bastante díspares, parece 
poder considerar-se como valor médio: 


Areias: v == 0,25 a 0,35 
Argilas: v == 0,45 a 0,50 


O módulo de Young, em terrenos arenosos, 
não depende da humidade nem da porosidade, 
variando pouco com a granulometria da areia. 
Apresentam-se valores médios para várias gra- 
nulometrias, no quadro 1: 


QUADRO 1 


Módulo de Young para terrenos arenosos de granulo- 
metrias diversas 


Granulometria E kg cm? 
(mm) (Módulo de Young) 
1,26 — 1,50 450 
1,00 — 1,25 520 
0,60 — 0,80 620 
0,35 — 0,60 480 
0,80 — 0,35 480 
0,20 — 0,30 620 


Em terrenos coesivos o módulo de Young de- 
pende considerávelmente das propriedades refe- 
ridas. Assim : 

Variação com o teor de humidade: pode ser ex- 
pressa pela relação 

E Es Eo (1 us as 
Wo? 
Eo — módulo de Young em solos de humidade 
w==0 

Wo — teor de humidade para o qual E = 0 

Variação com o teor de areia: é praticamente 
nula. O valor wo, porém depende deste teor. 
Assim, para argilas com menos de 50 / de areia, 
Wo varia na razão inversa da percentagem de 
areia. Por exemplo, se para uma argila pura, 
wo = 29"/o, para argila com 20º/, de areia, wo == 
= 20"/; para argila com 30 a 40'/y de areia, 
wo==20"/ e finalmente para argila com 50 a 
60º/y de areia, wo = 14 a 15º/. 

Para argilas com mais de 60º/o de areia, o mó- 
dulo de Young não varia com wo. 

Variação com a porosidade: Podemos considerar 
uma relação linear muito sensível. Por exemplo, 
numa determinada argila a variação da porosi- 
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dade de 0,32 para 0,93 conduziu a uma varia- 
ção do módulo de Young de 1700 para 120 
kg/cm*. (1) | 

Apresenta-se no quadro 2 valores do módulo 
de Young obtidos para vários tipos de solos. 
Note-se que para a resolução de um problema 
particular é indispensável a respectiva determi- 
nação por via experimental. 


QUADRO 2 
Módulos de Young de alguns solos 


Solo E, kg/cm? 

Argila plástica, com areia e silte 10 
Argila castanha, saturada, com 

silte e areia 440 
Argila densa com silte e alguma 

areia 2,950 
Areia de humidade média | 540 
Areia cinzenta com gravilha | o4U 
Areia fina, saturada 850 
Areia corrente 830 
Loess, areolas 1000 — 1300 


1.3 — Coeficiente de compressão elástica uni- 
forme, Cu 


Submetido a uma carga aplicada através de 
uma placa rígida, o solo sofre um assentamento. 
Retirada a carga verifica-se que a restituição é 
parcial. O assentamento é pois constituído por 
uma parte elástica e por outra residual. 

A relação entre a pressão e o assentamento 
total exprime-se através de gráficos com o aspecto 
da tig. 3. Nota-se um troço linear seguido dum 


e 


Pressoo 


Assentamento, mm 


Fig. 3 — Ensaio de carga sobre areia, com uma placa 
rígida de 1,0 m?, (!) 
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Pressão 


troço não linear em que o assentamento aumenta 
rapidamente. 

O limite do 1.º troço, linear, normalmente é 
tomado como critério para a capacidade de carga, 
estática, do terreno. 

Verifica-se que este limite depende apenas das 
propriedades do solo, sendo independente da área 
da placa de carga. 

Abaixo deste limite o solo pode considerar-se 
portanto como um sólido de deformação linear. 


Assentamento 


Fig. 4 — Ensaio de carga e descarga num solo. 


O assentamento elástico, de grande importân- 
cia no cálculo de vibrações nos solos, varia linear- 
mente com a pressão máxima em cada ciclo, 
segundo a relação 


o Cu de (3) 
sendo pz — pressão aplicada 
de — assentamento elástico 
Cu — coeficiente de compressão elástica uniforme 


O coeficiente de compressão elástica uniforme 
depende não só das propriedades do solo, como 
da área da placa de carga. 

Para uma placa rígida circular de área A, tem 
o valor: 


E 1 
Ê 1-» VA 
sendo pois 
E de 
= 1,13 + SRS (5 
E 1-) VA 


Na prática o valor Cu difere muito pouco do 
coeficiente de proporcionalidade entre a pressão 
e o assentamento total, devido às tensões de 
segurança sob os blocos de fundação terem sem- 
pre valores reduzidos. 


84 


A pressão varia com a distância 2, ao centro 
da placa, segundo a relação 


E Fa 


A eq 6 
1—» nVRê—p? (6) 


Pp — 


sendo R o raio da placa. 
Se a placa de carga for absolutamente flexível 
a pressão é uniforme e o assentamento variável. 


Designando por 4 o assentamento médio, de- 
monstra-se que para uma placa de secção rectan- 
gular de lados 2a e 2b, o valor do coeficiente de 
compressão elásiica uniforme é dado por 


Cs. € 
Cu = (7 
u E ) 
sendo 
C= Ja 
1—y 
e Cs função da relação «x = E , Conforme se 


indica no quadro seguinte : 


QUADRO 3 
Vulores de Co e C'. 


a | Ca C'. 
o —— 
1 | 1,06 1,US 
1,5 | 1,07 1,09 
| 1,09 IO 
2 1,18 | 1,15 
3 122 | 1,24 
10 | 1,41 1,42 


Os valores C', referem-se a placas rígidas. Na 
prática diferem pouco de C.. 

Apresenta-se no quadro 4 valores indicados 
pelo Regulamento Russo de Cálculo e Cons- 
trução de Fundações de Máquinas para o coefi- 
ciente C,, para áreas de fundação de 10 mº, 

É de notar que a variação real de Cu com a área 
de fundação afasta-se um pouco da teórica, pelo 
que é perigoso extrapolar valores de C, a partir 
de ensaios sobre placas de dimensões muito pe- 
quenas (como é o caso dos ensaios edométricos). 
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QUADRO 4 


Valores do coeficiente de compressão elástica uniforme, Cu, para áreas de 10 mº 


Capacidade 
Categoria de carga | Coeficiente 
do Descrição do solo estática Cu 
solo Kelem? Kejem? 
Solos fracos: argilas, plásticas, com silte 
ou areia; areias com silte ou argila; À 
solos das categorias II e III com silte | Até 1,6 | Até 5 
orgânico. 
H Solos médios : argilas, argilas com silte LBlBD] pos 
e areia, perto do estado plástico ; areia. 
Solos fortes: argilas sem ou com silte 
HI ou areia, compactas; gravilha; areia | 3,9 —5 5 — 10 
com gravilha; loess; areolas. 
IV Rochas pd, > 10 
1.4 — Coeficiente de compressão elástica não sendo 
uniforme, Cu C= 1 
J— vê 


Seja uma placa rígida assente sobre o solo 
sujeita a um momento M que a faz rodar em 
torno do eixo principal. Em cada ponto da placa 
desenvolve-se uma reacção vertical, elementar, 
devida à resistência do solo, dada por: 


dR = C; Lo dA (8) 
sendo : 


dA — elemento de superfície da placa 

7 — ângulo de rotação 

L = distância de dA ao eixo de rotação 
Cy — coeficiente de compressão elástica não uniforme 
dR — reacção elementar 


O momento de dR em relação ao eixo será dM 
dM = C; Lº 92 dA 
e o momento total M 


M=C,p/L'dA=C.9.] 
A 


sendo I o momento de inércia da base da placa 
em relação ao eixo. 

Para uma placa rectangular de lados 2a e 2b, 
demonstra-se ser. 


Cy —— = (10) 
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e K função da relação « == A conforme se in- 


dica no quadro 5. No mesmo quadro apresenta-se 
a relação 


Cs Ko 
- (a) = (11) 
Cu Cs 
QUADRO 5 
Cu 


Valores de Ky e —— 


n 


É Et 
LO | 1,984 | 1,87 
1,5 | 2,254 | 2411 
2.0 | 2,510 
3 |2,95| 2,68 
5 | 3,700! 8,04 

IO | 4,981 


Verifica-se que C, é sempre inferior a C7 po- 


4 * - Co 
dendo tomar-se em média a relação —— =2 
u 
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QUADRO 6 


Valores dos coeficiente Cu e C; para diferentes tipos de solo e áreas de fundações 


| | Área da placa | Cu | Cy | Co 
| E sia | de carga (m)* | kg em* | kg/em Cu 
nte TS E a i a => mm 
| mm 
Argila castanha, saturada, 2,0 4,40 12,0 2,13 
com silte e alguma areia 4,0 2,50 4,0 1,60 
8,0 2,05 3,0 1,46 
Argila cinzenta, saturada, 0,5 3,5 3,50 1,02 
branda, com silte e alguma 1,0 2,52 | 3,61 1,d4 
areia 1,5 2,11 | 3,79 1,80 
Loess com humidade média 0,81 14,2 25,0 1,76 
1,4 10,8 17,6 1,63 
2,0 10,2 15,5 1,51 
4,0 8,0 12,9 1,61 
Areia cinzenta saturada, 4,0 T,o 14,5 1,92 
densa e fina | 8,0 2,6 9,5 IT 
15,0 4,0 9,2 2,80 


Apresenta-se no quadro 6 valores experimen- 
tais dos coeficientes referidos, para diferentes 
tipos de solo (!) 


1.5 — Coeficiente de corte elástico 
uniforme, C: 


Seja uma placa rígida assente sobre o solo, 
submetida à acção de uma força horizontal. A 
tensão média de corte aplicada, 7 e o respectivo 
deslocamento elástico transversal “e estão rela- 
cionados por 

t=C: Ae (12) 
sendo : 


C, = corficiente de corte elástico uniforme 


Para uma placa rectangular de lados 2a e 2b, 
demonstra-se ser 


K: € 
C = — (13) 
“VA 
sendo KZ Função da relação « = à e do coe- 


ticiente de Poisson, », conforme indica o quadro 7. 
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Demonstra-se ser 


Cu 


(14) 
C, K,; 


Como o coeficiente de Poisson » é maior em 
solos coesivos que em areias, O valor do coefi- 
ciente C. é maior em solos arenosos que em 
argilas. 

Verifica-se que C; aumenta proporcionalmente 
à pressão normal sobre o solo, e aumenta com a 
duração da carga aplicada. 

Verifica-se ainda que para áreas superiores a 
cerca de 12 m”, C; mantém-se constante com a 
área. 

Dados que os valores teóricos deduzidos para 
a relação Cu /C, não entram em linha de conta 
com a influência da pressão normal sobre o valor 
de Cu, aconselha-se tomar para esta relação os 
valores experimentais obtidos para casos corren- 
tes, de valor médio. 


Cu 
C, 


= 2 (14) 
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QUADRO 7 


Valores de KT em função de « e »v 


o Sa | 0,5 1,0 | 1,5 
041 1,040 | 1,000 | 1,010 
0,2 0,990 | 0.988 | 0,942 
0,3 | 0,926 | 0,868 | 0,864 
0,4 | 0,844 | 0,792 | 0,770 
0,5 | OITO | O,704 | 0,692 


1.6 — Coeficiente de corte elástico não uni- 
forme, Cy. 


Se uma fundação rígida assente sobre o solo 
for actuada por um momento Mz em relação a 
um eixo vertical, rodará em torno deste eixo 
de um ângulo 4 dado por: 

J. — momento polar de inércia da base da 
tundação 
C: — coeficiente de corte elástico não uniforme 


Verifica-se que o valor de C: é ligeiramente 
superior ao de C-, apresentando-se no quadro 8 
valores dos 2 coeficientes. 


QUADRO 8 
Valores de Cy e €: 
: Cx, Cu; SP — 
Tipo de solo A, m? ki/cnit eae? Cy/C- | 
Argila cinzentasa- | 10) 1,20 | 2,92 | 2,14 
turada, comssilte| 1.0| 1,64| 2,55| 1,42 
e alguma areia 0,5 | 1,88 | 2,08| 1,01 


DD —— ——— 


1,0 2,04 | 8,90 | 1,54 
Areia cinzenta fina 


4,0 | 2,20 | 285 | 1,30 
saturada 
15,0 | 1,90 | 2,08 | 1,16 
Aconselha-se a tomar em média 
| = 1,5 (16) 
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20 | 3,0 | 5,0 | 10,0 
1,020 | 1,050 | 1,150 | 1,250 
0,945 | 0,975 | 1,050 | 1,160 
o,870 | 0,906 | 0,950 | 1,040 
0,784 | 0,806 | 0,850 | 0,940 
0,586 | O,700 | 0,752 | 0,940 


1.7 — Resistência elástica de fundações 
sobre estacas 


1.7.1 — Acção de cargas verticais axiais 


O emprego de estacas em fundações sujeitas 
a choques e vibrações pode ser necessário em 
quatro casos : 


— À tensão total sobre o solo, estática e diná- 
mica, excede a capacidade de carga do solo; 

— é necessário aumentar a frequência própria 
das vibrações da fundação; 

— é necessário diminuir a amplitude de vibra- 
ções livres ou forçadas ; 

— é necessário diminuir o assentamento resi- 
dual, dinâmico, da fundação. 


O segundo e terceiro casos referem-se princi- 
palmente a fundações da máquinas com solicita- 
ções dinâmicas; os restantes são correntes no 
estudo de fundações de estruturas sujeitas à in- 
fluência dinâmica de ondas provenientes de ou- 
tras máquinas, propagadas através do solo. 

No primeiro caso, aplicam-se os métodos clás- 
sicos de determinação da capacidade de carga 
das estacas. No segundo e terceiro, surgem algu- 
mas particularidades, que se referem em seguida: 

Assim, se uma estaca é submetida a uma série 
de cargas e descargas, verifica-se que o assenta- 
mento elástico ô é proporcional à amplitude da 
carga P, sendo 

P=C; Ó (17) 
Cs — coeficiente de resistência elástica da estaca 


É fácil determinar por via experimental o va- 
lor de C5. Assim, a frequência das vibrações do 
solo produzidas pela cravação da estaca são 
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+. “ 1 A “ A . 
iguais à frequência própria da estaca, sendo 


T=271/ ê ss (18) 
V Cs 

T — período das vibrações 

m — massa da estaca 


Verifica-se que durante a cravação, a profun- 
didade da estaca não tem efeito sensível no 
valor de T, ou seja, a resistência elástica exer- 
ce-se principalmente através da base da estaca. 

Com o tempo, aumenta a coesão entre o solo e 
a superfície lateral da estaca, aumentando consi- 
deravelmente a capacidade de resistência elástica. 

O coeficiente de resistência elástica de um 
grupo de n estacas será 


Kz 


n Cx (18) 


sendo 


C;=v Cs 


“ — coeficiente <" 1, variável com a distância 
entre estacas, 

Sendo d o diâmetro das estacas e t a distância 
entre elas, é corrente ser: 


Apresenta-se no quadro 9 valores de €,, obti- 
dos experimentalmente por Barkan (1). 


1.7.2 — Acção de cargas verticais excêntricas 


Sejá uma fundação, assente sobre um conjunto 
de n estacas, actuada por um momento M, ro- 
dando consequentemente de um ângulo 7 em 
torno dum eixo principal horizontal da base. 
Sendo Li a distância da estaca i a esse eixo, será: 


n 
M = C&' 9º 2 L;* 
i=] 
ou: M = ks. ? (19) 
k,=C: à L; 


k. — coeficiente de compressão elástica não 
uniforme do grupo de estacas. 


Verifica-se que o valor de €; é muito seme- 
lhante ao de C5. 


1.7.3 — Acção de cargas horizontais 


Actuadas pos forças horizontais aplicadas no 


= 6; u = 0,65 
topo, as estacas comportam-se como colunas en- 
QUADRO 9 
Valores de C;; 
Área da Nºd Ga Comp das p A | Es 
e 
Tipo de máquina Tipo de solo fundação | E seria E asno ni 
m? | m m | Kg nt 
| 
Ra — Es apena asd 
Vindas Argilas brandas into 15 ass 0.55 3.0 0.6 0,1><10º 
com silte e pouca areia ps 
Ser A . . 4 . 4 
ra mecânica | Argila branca, Plástica, 17.4 10 | 0,50 65 08 1.3>< 10! 
de madeira com areia média | 
Areia saturadas com silte | | 
» oruânico e argila plás-| 24,5 30 ITU | 6,0) 0,8 3,4><105 
tica azul 1 | | | 
it: Bi reste ara | ao la! am | má lagoa 
» , á 24,8 20 0,80 6,0 0,8 TA 
turada, com silte 
. R 12,0 | 35 |0,40| 30 | 08 |5,3>x<10! 
a ESSE | AO dn, a | Becas 
» | » 18,0 go | 0,50 3,0 | 0,8 |87 S><10* 
Vibrador Loess 4,0 | y U,7D 6,U | 0,8 1, 1>< 10% 
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ANTIGAMENTE : 


HOJE : 


Com 
imagem 
direita 

na luneta 


Trabalho cómodo e simples 
sem erros e sem dúvidas, 
sinais correctos, não inver- 
tidos, claros e rápidos. 


O pequeno nível WILD N 10 com a 
designação «E» tem imagem direita do 
campo de observação na luneta. À mesma 
tem luminosidade igual ao modelo con- 
vencional e o mesmo tamanho. 


Representantes exclusivos 
Á 
WILD PORTUGAL, LDA. úissoaz — reLerone cs ni27 
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castradas a uma profundidade média de 1 a 1,50 
m. Verifica-se pois que não oferece vantagem 
elevar muito a altura da estaca. 

Assim, o factor preponderante na resistência 
da estaca a forças horizontais é a sua rigidez de 


flexão. 
Não é portanto aconselhável empregar estacas 


para aumentar a frequência própria de funda- 
ções, para vibrações horizontais. 


2 — ASSENTAMENTOS RESIDUAIS DEVIDOS 
A VIBRAÇÕES 


21— Introdução 


Verifica-se que fundações sujeitas a choques e 
vibrações são em certos casos susceptíveis de so- 
frer assentamentos consideráveis, uniformes ou 
não, que podem comprometer seriamente a se- 
gurança das estruturas. 

Especialmente no caso de solos não coesivos, 
geralmente com grande capacidade de carga para 
solicitações estáticas, a acção de solicitações em 
regime dinâmico leva a uma alteração conside- 
rável nas propriedades do solo, nomeadamente 
nas forças de atrito, de coesão, nas proprieda- 
des hidrodinâmicas, nos módulos de elasticidade 
e plasticidade. Esta alteração pode ser tal que o 
solo adquira características semelhantes às de um 
líquido viscoso. 

Seguem-se algumas indicações, de carácter 
essencialmente qualitativo, sobre a influência das 
vibrações nas propriedades enunciadas. 


2.2 — Influência das vibrações nas caracteris- 
ticas de amortecimento dos solos. 


22.1 — Atrito interno e coesão durante as vibrações 


Verifica-se que o coeficiente de atrito interno 
varia com a energia cinética das vibrações, po- 
dendo ser reduzido para 60 a 70º/y do seu va- 
lor inicial. Este parâmetro é de grande impor- 
tância na resistência ao corte. Aquela redução é 
menor em solos húmidos que em solos secos. Em 
areias, é mínima para um teor de humidade de 
12 a 13º/. 

O coeficiente de atrito interno diminui qando 
aumenta a amplitude ou a trequência das vibrações. 
Enquanto que a variação com a amplitude é re- 
gular, existe para cada solo uma banda de fre- 
quências para a qual não há efeito sensível sobre 
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o atrito interno; com o aumento da frequência 
para além do valor limite daquela banda, o atrito 
diminui bruscamente. 


2.2.2 — Amortecimento viscoso 


Considera-se que o solo se comporta como um 
sólido elasto-plástico. A análise dos respectivos 
ciclos de histerese sob a acção de cargas está- 
ticas e dinâmicas fornece indicações sobre o 
amortecimento. 

Sendo AW a área de um ciclo de histerese, 
proporcional à energia absorvida, e W o traba- 
lho total realizado até à máxima deformação, de- 
fine-se 


4 — coeficiente de absorção de energia 
Em sistemas com um grau de liberdade, em 
vibrações livres, é 


4=1—-erunT 


n — amortecimento viscoso relativo 
T— período das vibrações. 


Verifica-se que o coeficiente > de um dado solo 
é relativamente estável. Assim, varia pouco com 
o teor de humidade. É menor para compressões 
que para esforços de corte. Nas areias, varia con- 
siderávelmente com a granulometria, conforme 
se indica no quadro 10. 


QUADRO 10 


Variação de b com a granulometria, em areias 


Granulometria 


mm 
0,10 — 0,25 0,64 
0,25 — 0,50 0,68 
1,00 — 2,09 0,79 


2.3 — Infiuência das vibrações na porosidade 
e propriedades hidráulicas dos solos 


A influência das vibrações sobre estas pro- 
priedades é de grande importância tanto sob o 
ponto de vista de cálculo e desenho de funda- 
ções como sob o do emprego de vibrações para 
melhorar as características dos solos. 
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Em solos não coesivos, verifica-se que o grau 
de compactação produzido por vibrações depende 
só da aceleração destas. 

Em todos os solos a relação entre o grau de 
compactação (traduzido pelo índice de vasios) e 
a aceleração das vibrações é do tipo indicado na 
fig, 5. 


Indice de vazios 


4 


Aceleração 


Fig. 5 — Relação entre o grau de compactação 
e a aceleração das vibrações 


Verifica-se que o teor de humidade tem grande 
influência no andamento da curva da fig. 5. Em 
areias, o teor de humidade que permite a compac- 
tação máxima é da ordem dos 17 º/e. 

A compactação de um solo sob a acção de vi- 
brações produzidas por uma máquina pode levar 
assim a assentamentos residuais consideráveis, 

Verifica-se que a permeabilidade de um solo 
aumenta considerávelmente por acção de vibra- 
ções, o que se aproveita na prática para fazer a 
estabilização química dos solos (ex. por injecção 
de silicatos). 


2.4 — Elementos da teoria de assentamentos 
residuais de fundações sob a acção de 
solicitações dinâmicas 


Ainda que isso involva em certos casos erros 
consideráveis, o modêlo correntemente empre- 
gado no estudo de assentamentos residuais é o 
do oscilador elasto-plástico com amortecimento 


a 
18 
14 
14 
q 
- 
Q.. 
P| 
Assentamento 


Fig. 6 — Característica de restituição dos solos 
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viscoso, com a característica de restituição bi- 
-Jinear do tipo que se indica na fig. 6. 

Designando por pst e pd a tensão estática e di- 
nâmica sobre o solo para que não se produzam 
assentamentos residuais deverá verificar-se a 
relação 

pst + pd < pl (21) 

sendo pl a pressão admissível 

Designando por Rst e Rd respectivamente a 
tensão estática e dinâmica limites do assenta- 
mento residual, verifica-se que a relação entre 
pst e pd é do tipo da curva da fig. 7. 


Roá 


Pd B Ps f 


Rst 
Fig. 7 — Redução da pressão admissível, em presença 
de vibrações 


Substituindo a curva pela recta a tracejado, a 
condição (21) expressa-se por: 


Pst e q Rst 
sendo (22) 
2 pd 
E 1.== = 
' Rd 


sendo q o coeficiente de redução da tensão admis- 
sível do solo devido à presença de vibrações. 
Portanto, enquanto a pressão sobre o solo for 
menor que o valor limite pl, a fundação terá 
apenas assentamento elástico. 
No caso de oscilações verticais, sendo 


O Ma Cu Az 


a pressão dinâmica correspondente à amplitude 
de oscilação Az, obtém-se, a partir da inequação 
(21) a condição para que não haja assentamento 
residual: 


(23) 


Na maioria dos casos, desde que as fundações 
estejam bem dimensionadas, as amplitudes são 
muito inferiores a este limite, pelo que não ocor- 
rem assentamentos residuais importantes. 
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A determinação deste tipo de assentamentos 
assume grande interesse prático no caso de fun- 
dações sujeitas a oscilação livre provocada por 
choques verticais ou pela propagação de ondas 
através do solo. A amplitude destas oscilações 
pode atingir valores consideráveis, dando origem 
a pressões dinâmicas de valor elevado. 

A fig. 8 representa o diagrama da elongação 
vertical (assentamento-tempo) de uma fundação 
sujeita a um choque 


Fig. 8 - Elongação de oscilador elasto-plástico com amor 
tecimento viscoso, sujeito a um choque 
Zst — assentamento estático AB — assentamento elástico 


BC = Zo — Zst — assentamento residual 


À equação que rege o movimento do sistema 
na fase elástica será: 
Zi + 2CZ + w,? Zi=g (24) 
sendo : 


Z1 — assentamento elástico 
C — amortecimento 

vz — frequência própria 

g — aceleração da gravidade 


Quando a pressão atingir valores superiores 
ao limite elástico a força de restituição passa a 
ser constante e o movimento do sistema em fase 
plástica é dado pela equação : 


A+aCh=s(1— 5) (25) 


Pst 
Zs — assentamento residual 


O tratamento analítico deste movimento é 
simples e permite, dadas as características diná- 
micas do sistema e as condições iniciais, determi- 


nar completamente as funções Z (t), Z(t) e Z (t). 
A velocidade inicial da fundação para a qual 

começa a ocorrer assentamento residual é dada 

por 

Pi — Pst / 8 


eae (26) 
V Pst Cu 


Vo =) 
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sendo : 
j=5(0);0=C/m 


e estando j e p relacionados como se indica no 
quadro 11. 
QUADRO 11 


Valores j (9) 


1,00 
0,3 1,47 
0,6 1,84 


Para valores de velocidade inicial superiores 
a Vo, O assentamento residual é aproximada- 
mente proporcional âquela velocidade. 

Torna-se possível portanto dimensionar uma 
fundação de modo a evitar assentamentos dinâ- 
micos residuais importantes. 

O assunto será desenvolvido noutro artigo 
desta série. 


2.5 — Métodos para limitar o assentamento 
residual dinâmico de uma fundação 


Como se viu, a tensão admissível num solo, 
desde que se verifiquem vibrações, deve ser 
reduzida. Esta redução varia com as amplitudes 
de vibração. Assim, nos casos mais desfavorá- 
veis, que são os correspondentes a fundações de 
máquinas de percussões, tipo martelos, o valor 
do coeficiente de redução deve ser 


q=0,4 a 0,5 


Nestes casos a amplitude atinge correntemente 
valores de 1 a 1,5 mm. Para máquinas alterna- 
tivas, de frequência baixa e amplitudes que não 
excedam 0,2 a 0,3 mm, já o valor de q pode 


ser 
q = 1,0 


Em máquinas rotativas de alta velocidade (tur- 
binas) as amplitudes são da ordem da dezena de 
microns. No entanto, como estas máquinas são 
muito sensíveis a vibrações deve tomar-se um 


valor de 
q=0,8 


No caso de vibrações propagadas através do 
solo, a fim de evitar assentamentos residuais não 


uniformes de valor considerável, em fundações 
de outras máquinas ou edifícios, as tensões de 
segurança devem variar na razão inversa da dis- 
tância da respectiva fundação ao foco de propa- 
gação das ondas. 

É aconselhável no caso de terrenos arenosos 
sujeitos a grandes assentamentos, proceder à sua 
estabilização. Em regra é suficiente fazê-la até à 
profundidade de 3 ou 4 m. 

Outro processo de reduzir assentamentos diná- 
micos é o emprego de estacas curtas. Recomen- 
da-se o emprego de estacas de madeira, pela sua 
economia e eficácia. 


3. ELEMENTOS DE TEORIA DAS VIBRA- 
ÇÕES 


3.1. Generaiidades 


A fundação da maior parte das máquinas de 
média e grande potência é constituida por um 
bloco muito rígido assente ou embebido no solo, 
O funcionamento da máquina comunica ao bloco 
de fundação solicitações, periódicas ou não, que 
originam diversos tipos de movimentos oscilató- 
rios. Assim, o bloco da fig. 9 admite 6 tipos de 
movimento, correspondentes aos respectivos 6 
graus de liberdade: 3 translações e 3 rotações. 


Fig. 9 


Modos de vibração de um bloco rigido 


O sistema corresponde portanto a um oscila- 
dor de 6 graus de liberdade, assente sobre um 
meio parcialmente elástico, o solo, que influencia 
a dinâmica do bloco devido às suas característi- 
cas de massa, elasticidade e propriedades de dis- 
sipação de energia, ou amortecimento. 

De facto, verifica-se que uma parte do solo 
adjacente ao bloco vibra em fase com ele, aumen- 
tando a massa total em movimento. 

A análise dinâmica do sistema de 6 graus de 
liberdade assente sobre um meio parcialmente 
elástico semi-indefinido não está ainda formulada. 
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Reissner (2) apresentou em 1936 a resolução do 
problema para oscilações verticais de um bloco 
de base circular assente sobre um meio semi-in- 
definido, perfeitamente elástico. Os seus traba- 
lhos foram continuados e desenvolvidos por in- 
vestigadores de vários Países, entre os quais 
Arnold, Bycroft, Warburton, Shekhter, (4),(8) 
conduzindo ao estudo da propagação de ondas 
de diversos tipos no interior do meio elástico. 

Verifica-se porém que salvo em casos de pro- 
ximidade de ressonância, é lícito assemelhar o 
solo a uma mola sem massa, com característica 
de restituição linear, visto que as soluções assim 
obtidas não se afastam muito das reais. A aná- 
lise dinâmica dos blocos de fundação torna-se 
portanto muito simples, desde que se admitam 
as seguintes hipóteses: 

— Considera-se uma relação linear entre assen- 
tamentos elásticos e as respectivas forças de res- 
tituição. Esta relação é expressa pelos coeficien- 
tes de compressão e corte elásticos, uniformes e 
não uniformes, definidos em 1.3, 1.4, 1.5 e 1.6; 

— despreza-se a inércia do terreno; 

— considera-se que a fundação, no caso de 
ser constituída por um bloco maciço, tem uma 
rigidez muito elevada em relação ao solo; 

— desprezam-se as reacções do solo sobre a 
superfície lateral da fundação, quando esta é 
enterrada. Estas reacções têm contudo influência 
sensível sobre a frequência própria e o amorte- 
cimento, em regime livre. Devido a isso, no es- 
tudo de fundações de máquinas de percussão 
considera-se um acréscimo no coeficiente de res- 
tituição do solo, na base da fundação. 


3.2 — Vibrações verticais 
3.2.1 — Vibração livre sem amortecimento 


Dentro das hipóteses referidas, o movimento 
de um bloco sujeito a vibração livre, assente 
sobre o solo, é dado por: 


Z"'+ E Z=0 (27) 
sendo : 
Z — elongação vertical 
ft, — frequência própria tal que 
9 Cu. A 
Ê; — E (28) 
m 


C, — coeficiente de compressão elástica uni- 
forme do solo 


TECNICA Nº 852 


A — área da fundação 
m — massa do bloco e da máquina 


A solução de (27) é, como se sabe 


Z=aAsenf,t + Bcos f,t 


sendo A e B constantes dependentes das condi- 
ções iniciais. No caso particular de um choque 
será : 

t=0; Z,=0; Z0= Vo; 


Z = - sen fz t (29) 
Z 
vo — velocidade inicial 
3.2.2 -- Vibração forçada sem amortecimento 
O movimento é regido pela equação 
" P(t 
pt 
m 
sendo P (t) — a solicitação imposta 
No caso de ser 
P(t) =P sen ot 
a solução forçada é dada por 
P 
z= Azsenot= ss sen ot 
m (£, — (3º) 
Pp 1 1 
A=—s. 3 = Àst.—sa 0 
nA ioga top 


Az — amplitude de oscilação forçada 
Ast — assentamento para uma carga estática 


de valor F. 
» 


Na prática, sempre que empregam dispositivos 
para isolamento vibratório de fundações, dá-se- 
-lhes a menor frequência própria possível, o que 
reduz a amplitude 4.. 

Em máquinas alternativas, com velocidade 
baixa, convém que a frequência própria das fun- 
dações seja muito elevada. Note-se que para con- 
seguir isto, os blocos de fundação devem ter uma base 
grande e altura mínima, o que contraria a convic- 
ção corrente de que «a melhor fundação é a 
mais pesada”. 
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3.2.3 — Vibração com amortecimento viscoso 


Traduzido o amortecimento do solo pelo coefi- 
ciente de atrito C, a equação de um bloco de 
fundação assente sobre o solo será dada por 


z+2Cz+fz=0 (31) 
cuja solução é, com se sabe 
mo SÉ 867" = 
t=t “(Aseny/ Bp— Ct + 
(32) 


+ BcosVf!— C! t) 
sendo A e B constantes dependentes das condi- 
ções iniciais. 

No caso corrente de ser 


EC < f 


a equação (32) pode escrever-se, sem erro apre- 
ciável sob a forma 


PR (A sen f, t + Bcos f, t) 


3.2.4 — Vibração forçada com amortecimento viscoso 


No caso de solicitação sinusoidal do tipo 


P (t) =Psen mt 
será 


z +2Cz t+ ftz=Psenwt (33) 


cuja solução forçada é 

As = — E 
my (But) + 4Cio? 

tg y=2Co/(E — a?) 


Zz== À; sen (mt —:y) 


(34) 


Como se sabe, o amortecimento tem grande 
influência no valor da amplitude em zonas pró- 
ximas da ressonância. 


325 — Influência da inércia do solo em vibrações 
forçadas verticais 


Verifica-se que a influência da inércia do solo 
se traduz por um efeito equivalente a um au- 
mento de massa da fundação e uma dissipação 
de energia dependente não só das características 
do solo (e especialmente do coeficiente de Pois- 
son, v), como também da massa e dimensões da 
fundação (1), 
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Para um solo caracterizado por um coeficiente 
de Poisson v, definindo o coeficiente 


ade ps 


(35) 
/ VA 


Y — peso específico do solo 
Pst — pressão estática 
A — área da fundação 


determinam-se a partir dos gráficos da fig. 10 
os coeficientes de amortecimento é e de aumento 
de massa f. 


0 5 o b 15 20 


Fig. 10 — Valores de £ e 


A partir destes valores, a amplitude 4. das 
vibrações do bloco, sujeito a uma solicitação 
sinusoidal P sen mt é dada por 


V [3,84/(1-v) —x?bB]º 4 15,35 xºbB:2/(1—2) 


A /384(1—28/6) 
pet b(1—v) E 


G — módulo de elasticidade transversal do solo 
to—VA/R 
A — área da fundação 


Verifica-se que valores pequenos de b corres- 
pondem a fundações de base grande e altura 
reduzida. Estas fundações são as mais influen- 
ciadas pela inércia do solo. 

Na prática corrente pode desprezar-se o efeito 
da inércia do solo, por não conduzir a grandes 
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alterações nos valores da amplitude e frequência 
própria de vibrações, para blocos assentes no 
solo. No caso de blocos embebidos, já o efeito 
da inércia do solo pode ser considerável. 

Desde que a frequência forçada se aproxime 
da frequência própria do sistema bloco-máquina 
ou seja, na proximidade da ressonância, é indis- 
pensável considerar a inércia do solo no cálculo 
da amplitude de oscilação. 

É de notar que embora se possa quase sem- 
pre dimensionar as fundações de forma a afas- 
tar a sua frequência própria da velocidade de 


fd 


1 em função de v e b 


regime da máquina, ocorrem frequentemente 
fenômenos de ressonância, se bem que transitó- 
rios, nos períodos de arranque e paragem das 
máquinas. 


3.3 — Vibrações angulares em torno de um 
eixo horizontal 


Seja um bloco de fundação assente sobre o 
solo, simétrico em relação a um plano vertical, 
oscilando em torno do eixo horizontal (0), con- 
forme se indica na fig. 11. 

Seja W o peso do bloco. 


Fig 11 — Bluco sujeito a oscilação angular 
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Sendo 7 o ângulo de rotação do bloco em re- 
lação à posição média, a equação que rege o mo- 
vimento livre será 


Wo9 +(CGI-WL)s4=0 (37) 


sendo: 


Wo -——-momento de inércia do bloco em rela- 
ção ao eixo de rotação 

C, — coeficiente de compressão elástica não 
uniforme do solo 

| -——- momento de inércia da base em relação 
ao eixo de rotação 

L — cota do centro de gravidade 


A solução da equação (35) é 


2=C sen (fy t + 90) 
ph ss E VM 


(38) 

Wo 
f» —frequência própria da oscilação angular 
C, vo — constantes dependentes das condições 


iniciais 


Se o bloco estiver sujeito a um momento da 
forma 


M sen ot 


a solução forçada é, como se sabe 
senwt (39) 


Como geralmente é WL << Cy, a frequên- 
cia própria angular reduz-se na prática a 
o 

Wo 


fo? (40) 


Para fundações de base rectangular de lados 


a e b, sendo a o lado perpendicular ao eixo de 
rotação será 


ba co ba cs 
+ RE 122 Wo 


A dimensão a tem influência preponderante 
no valor da frequência própria. 

O assentamento máximo, na extremidade da 
fundação será 


add q Ma 


2 Na 2Wo (ES — 9º) 
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Analizou-se o caso de vibração sem amorte- 
cimento. É muito fácil introduzir o amorteci- 
mento viscoso nas equações. 

As oscilações angulares ocorrem principalmente 
em fundações altas de máquinas com componen- 
tes de forças horizontais ou de momentos não 
equilibradas. 


3.4 — Vibrações de corte puro 


« 


Se a resistência do solo à compressão for 
muito superior à resistência ao corte, ocorrem 
sobretudo movimentos horizontais da fundação, 
se existirem solicitações horizontais. 

Se o bloco for actuado por uma força hori- 
zontal 

P sen mt 


a equação de oscilação horizontal será 


lr 


x + fi x= 


E. sen mt (41) 
m 


x — elongação horizontal 


Ex — frequência própria da oscilação horizontal 


sendo 


(42) 


C, — coeficiente de corte elástico uniforme 
A — área da fundação 
m — massa do bloco 


A solução forçada de (39) é dada por 


-—— senwt 


A 43 
X m (4 — q?) À ) 


3 5-— Vibrações angulares em torno de um 
eixo vertical 


Se a fundação for actuada por um momento 
M senwt em relação a um eixo vertical, sendo 


W, — momento de inércia da fundação em re- 
lação ao eixo vertical 

J. —momento de inércia polar da base da 
fundação 

v — ângulo de rotação 

Cy — coeficiente de corte elástico não uniforme 
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será 


W,4 + Ci), y-=Msenot (44) 
sendo a solução forçada dada por 


— M 
W, (E, — (3º) 


É; =CiJ,/W; (45) 


| === 


sent 


3.6 — Vibrações compostas de deslocamento 
vertical, rotação e deslizamento 


Nos parágrafos anteriores apresentaram-se 
movimentos com 1 grau de liberdade. Estes mo- 
vimentos só são, em regra, possíveis, desde que 
o solo seja caracterizado por uma rigidez infi- 
nita em relação à compressão ou ao corte. Não 
sendo esse o caso real, a oscilação das funda- 
ções tem características mais complexas, corres- 
pondentes a movimentos com vários graus de li- 
berdade. 

Consideremos então um bloco simétrico em 
relação a um plano vertical, assente sobre o 
solo, actuado por uma força P (t), aplicada no 
centro da gavidade, G, do bloco e um momento 
M (t), existentes no plano vertical de simetria, 
conforme se indica na fig. 12. 


Fig. 12 — Bloco sujeito a oscilação composta 


Demonstra-se facilmente que o movimento 
deste sistema é dado pelas equações seguintes: 
mz +C Az=P,(t) 

mx + C Ax—-C, AL>2=Px(t) (46) 
Mn2C,ALx + (CI-WL + CAL)2=M (t) 


sendo ; 
m — —massa do bloco 
Z — elongação vertical 
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Ci —coeficiente de compressão elástica uni- 
forme do solo 

A — área da fundação 

P, (t) — componente vertical da força aplicada 

x — elongação horizontal segundo x 

C,  —coeficiente de corte elástico uniforme do 
solo 

L — cota do centro de gravidade do bloco 

o — rotação do bloco em torno do eixo hori- 


zontal que passa por G 
Px (t) — componente horizontal, segundo xr, da 
força aplicada 


Mm — momento de inércia do bloco em relação 
ao eixo horizontal que passa por G 

C. —coeficiente de compressão elástica não 
uniforme do solo 

I — momento de inércia da base em relação 


ao seu eixo horizontal principal 
W  —peso do bloco 
M (t) — momento aplicado. 


A 1.2 equação do sistema é independente das 
restantes. Portanto, se a solicitação não tiver 


+ 


componente vertical, é sempre z(t)=0. Sea 
única solicitação actuante for vertical, será 
x (t) =» (t) = 0. 

Analisada a 1.2 equação nos parágrafos an- 
teriores, consideremos o sistema definido pela 
24 e 3.4 equações. 


361 - Vibração livre 
Neste caso, é 
mx +CAx—-C,AL2=0 (47) 
Mas —C, ALI +(CGI-WL+C,AL)9=0 


o sistema tem soluções particulares da forma 


x = Aa sen (ft + a) 
9 = Ba sen (Ft + 2) 


Eliminando x 7 no sistema (45) obtem-se a equa- 
ção das frequências 


E — Cl — WI L C, À L*m + e E 4 
N Mm m Mm : 
pGI-WL GA, (ug 
Mm m 
Designando por Mo 
Vice St EE aj 
m 
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Existe uma borracha para o revestimento destes rolos que devem 
manter-se em pleno funcionamento a temperaturas até 180º €.? 


Também este problema foi resolvido com *Silopren! 

Os fabricantes de rolos recorrem sempre ao Silopren para 
revestir rolos que, em condições normais de funcionamento, 
tenham de ter grande resistência ao calor, e, em presença de 
materiais pegajosos, produzam um bom efeito antiaderente. 
Novos tipos especiais oferecem consideráveis vantagens na 
laboração, o que permite uma fabricação económica: podem ser 
vulcanizados tal como a borracha natural; não é necessário um 
aquecimento posterior. Dai resulta, portanto, a possibilidade de 
fabricar rolos de quaisquer dimensões. É possível conseguir uma 
boa aderência do revestimento ao núcleo de metal, mesmo sem 
empregar agentes de aderência, utilizando o Silopren SVA Batch 
||-SH, que tem aderência própria e que durante a vulcanização 
se une fortemente a muitos metais. Além de apresentar muitas 


outras propriedades vantajosas, tais revestimentos de rolos. 


distinguem-se pelo facto de continuarem em bom funcionamento 
mesmo a temperaturas até 180º C. Durante curtos espaços de 
tempo resistem até mesmo a temperaturas muito mais elevadas. 
Além disso são fisiolôgicamente inertes e produzem efeito repe- 
lente em relação a outros altopolímeros. 


Por isso preferem-se rolos com revestimento de Silopren para 
a termocolagem de artigos de plástico para os sectores de 
géneros alimentícios e de cosméticos; e ainda para o revesti- 
mento de papeis, metais e téxteis com plásticos (por exemplo 
com polietileno), em que se tem de contar com temperaturas 
elevadas. Na cerâmica e na indústria vidreira utilizam-se rolos 
de Silopren como rolos de transporte e de guia, assim como 
rolos para o arrefecimento de material quente. 

Se quizer saber mais pormenores sobre o Silopren, escreva, por 
favor, aos nossos representantes. Eles terão muito prazer em lhe 
mandar material de informação muito interessante. 


S.A.R.L. QUIMICOR, REPRESENTANTE EM PORTUGAL - 
LISBOA: RUA SOCIEDADE FARMACÊUTICA, 3, TEL. 42194 - 
PORTO: RUA DE SANTOS POUSADA 441. TEL. 54141 
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o momento de inércia do bloco em relação ao 
eixo horizontal que passa pelo centro da base e 
considerando os valores 


das frequências próprias para oscilação e corte 
puro, a equação (48) transforma-se em 


. 
É; . És 
-— t 
/ 7 


Esta equação tem 2 soluções positivas, fi e fa, 
correspondentes às frequências próprias do sis- 
tema, verificando-se a relação 


fa < fo, fx < fu 


No caso geral, o movimento do bloco é a so- 
breposição de oscilações de frequência fi e fs, 
sendo a forma do movimento definida pelo va- 
lor 


(50) 


Se a fundação vibra apenas segundo a fre- 
quência mais baixa, fz, A, e B, são do mesmo 
sinal, correspondendo o movimento a uma rota- 
ção em turno de O (fig. 13). 

Se a fundação vibra só segundo a frequência 
mais elevada, fi, Aa e Ba são de sinais contrá- 
rios e o movimento é do tipo indicado na fig. 14. 

4 e pa serão definidos pela equação (50), sub- 
stituindo f respectivamente por fi e fi. 


i 
| 
Fig. 13 — Movimento segundo a frequência 


própria inferior 
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Note-se que é 


(51) 


No caso geral, o movimento do bloco corres- 
ponde pois a uma oscilação em torno de um 
ponto que dista p do seu centro de gravidade. 


3.6 2 — Vibração forçada 


Seja o caso do sistema actuado por uma força 
horizontal Psenwt. As equações (44) tomam a 
forma 


mx + C;Ax— C-AL72=P sen qt (52) 
M wu? — C:ALx + (CH-WL +CG AL)9=0 
sendo as soluções forçadas dadas por 


x == As sen vt = 


— CoI-WL + CAL? — Mo? 


P sen ot 
A (192) 


C:AL 
2 P sen qt (53) 


v = Aysen wt = 


A (52) = m Ma (£2 — vs?) (£2— 62) 


Se o sistema for actuado por um momento 
M sen wt, será 
C,AL 
Ax —— : 
A (mm?) 


54 
CA—-mo 


Á (12) 


Fig. 14 -- Movimento segundo a frequência 
própria superior 
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Existem 2 possibilidades de ressonância. À 
cada frequência w corresponde um tipo de mo- 
vimento constituído por uma oscilação em torno 
de um eixo horizontal afastado da distância p do 
centro de gravidade do bloco, sendo 


— Ax 


== (55) 
A; 


y 
Note-se que: para A;=0, o movimento é só 
de translacção 
para Ax==0, o movimento é de 
rotação em torno do centro de 
gravidade do bloco. 


3.6 3 — Excentricidade do centro de gravidade em rela- 
ção ao centro da base. 


Em todos os casos anteriormente analisados, 
o centro de gravidade do conjunto bloco-má- 
quina e o centro da base estavam na mesma ver- 
tical. Se a distribuição de massas for excêntrica 
em relação à base, o que pode suceder no caso 
de 2 máquinas acopladas no mesmo veio, ou de 
irregularidades de forma na fundação, há altera- 
ções sensíveis nas características de oscilação. 

Seja o caso mais simples em que o centro de 
gravidade e o centro da base ficam no mesmo 
plano vertical principal, mas não na mesma ver- 
tical. Analisemos a oscilação no plano principal. 

Considerando o sistema de eixos coordenados 
indicados na fig. 11, demonstra-se que a oscila- 
ção livre do sistema é dada pelas equações 


mx +C,Ax- CALv= 


mz + CG Az-C Asp=0 (56) 


Mn9— C,ALx + (CoI— WL+ 
+CGAS+CAI)9—CA:Z=D 


adoptando a rotação indicada e sendo 


: — excentricidade do centro de gravidade em 
relação ao centro da base. 


Ao contrário do que sucede quando não há ex- 
centricidade (46), neste caso existe dependência 
entre as 3 equações (56). O sistema admite 3 
frequências próprias de vibração, e qualquer so- 
licitação dá origem a um movimento constituído 
pela sobreposição das sinusoides: 


x == Ás sen fit 
z=B,sen f: (t— 1) 
y = Ca sen f; (t — 23) 


Designando por 
C, A = CC; Cy A = C, 


ClI-WL+(Cº+ CGL) A=C 


Obtém-se a equação fundamental das frequên- 
cias próprias : 


(Cum) (...)Ge=mf) EC >-Ma FP) = 
— (Ck—-mf)Cié—[(C;—mf)CiLº=o 


cujas raízes fi, fz e f; são as frequências pró- 
prias do sistema. 

Devido à excentricidade, as 2 frequências in- 
feriores são menores que as frequências inferio- 
res no caso de excentricidade nula; a frequência 
superior é maior que a frequência própria má- 
xima no caso de excentricidade nula. 

Na prática, desde que a excentricidade não ex- 
ceda 5º) do comprimento do lado menor da 
fundação, o seu efeito pode desprezar-se. 


3.7 — Comparação entre resultados analíticos 
e experimentais 


371 — Vibrações verticais 


D. D. Barkan, A. Mikhalchuk e N. Smolikov 
realizaram uma série de ensaios estáticos sobre 
solos, para a determinação do coeficiente de 
compressão elástica uniforme. Sujeitando-os em 
seguida a vibrações verticais transmitidas por 
fundações de áreas variáveis, verificaram que a 
frequência própria da fundação dado pelo cál- 
culo não diferia, em média, mais de 10 º/, da 
frequência própria medida. 

Verificaram também que a influência da inér- 
cia do solo sobre a vibração era muito reduzida, 
para uma série de fundações correntes a que 
correspondia um valor de b 20 (3.2.5). 


3.7.2 — Coeficiente de amortecimento 


O valor do coeficiente de amortecimento dum 
solo está relacionado com a amplitude, 4., das 
vibrações verticais dum bloco de fundação, em 
ressonância, pela expressão 


E 
Nip e ae 
2m é f, 
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sendo Psenft a força aplicada e m a massa da 
fundação. 

É possível pois, medindo a amplitude das vi- 
brações em ressonância, determinar o valor $. 

Apresenta-se no quadro 12, os valores de b e £ 
determinados por Barkan para fundações de 
áreas diversas, assentes sobre diferentes tipos 
de solos. 


rar erros de 10 a 20 º/ nos valores de cálculo das 
amplitudes. 


3.7 3 — Vibração de fundações sobre estacas 


Barkan levou a efeito uma série de determina- 
ções do coeficiente Kz em fundações sobre esta- 
cas, calculando em seguida a frequência própria 


QUADRO 12 


Valorees de b e £ para vários tipos de solos 


Tipo de solo 


Argila castanha saturada de água, 
com silte e alguma areia 


Argila cinzenta branca saturada 
de água com areia e silte orgânico 


Areia cinzenta, fina e densa, satu- 
rada de água 


Área da Peso da 

q Ga rá b E 
ai Ton. | 
2 16,3 16,5 | 0,145 
4 18,8 6,7 | 0,138 
8 30,0 40 | 0,88 
0,5 325 | 260 | 0,071 
| Do 14,5 0,058 
1,5 6,72 9,0 | 0,05] 

| à | 676! 190 | 0,132 

] 6,76 20 0,190 
l 6,70 20 0,175 
1 16,4 55 0,083 


O coeficiente de amortecimento em argilas 
brandas cinzentas é maior que em argilas casta- 
nhas. E é maior em areias que em argilas cas- 
tanhas. 

Verifica-se que os valores do quadro 12 são 
sensivelmente inferiores aos valores analíticos, 
dados pela fig. 10. Afigura-se que ós valores da 
fig. 10 correspondem a fundações embebidas no 
terreno, enquanto que os valores do quadro 11 
se referem a fundações assentes sobre o terreno. 

Em média, a relação entre o amortecimento 
naquele e neste caso é de 3 a 3,5. A relação en- 
tre o amortecimento de fundações parcialmente 
submersas em águas e simplesmente assentes é 
de, em média 1,5 a 2. 

Por outro lado o coeficiente 4 aumenta com a 
amplitude das oscilações o que não é tomado em 
conta nesta análise. 

No caso de vibrações afastadas da ressonân- 
cia, desde que os coeficientes de compressão es- 
tejam correctamente determinados, são de espe- 
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da fundação para vibrações verticais. Compa- 
rando os resultados obtidos com ós determinados 
experimentalmente, verificou que a frequên- 
cia própria real é inferior à dada pelo cálculo, cor- 
respondendo a um oscilador de massa superior 
em cerca de 30º, à massa da fundação. 

Este acréscimo corresponde à massa do terreno 
que oscila em conjunto com a fundação. 
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ERRATA 


NUMERAÇÃO DOS VOLUMES DA TÉCNICA 


A partir do número 331 a numeração dos volumes da «Técnica» saíu errada (no frontispício 


da Revista e no índice dos volumes). 
Correcção a fazer : 


Números em vez de deve ler-se 

331 — 340 Volume XXV Volume XXVI 

341 — 350 Volume XXVI Volume XXVII 
351 Volume XXVII Volume XXVIII 


Do facto pedimos desculpa aos nossos leitores. 
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AQUECIMENTO DE 
IMERSOS 


(Conclusão) 


4 — CALOR PERDIDO 


A degradação de energia electromagnética em 
calor num transformador provoca o seu aqueci- 
mento. Estabelece-se assim entre o transforma- 
dor e o meio ambiente uma diferença de tempe- 
raturas que provoca a transferência desse calor 
ao referido meio, pelos bem conhecidos meca- 
nismos de condução, convecção e radiação. 

A condução é responsável pela transferência 
duma parte do calor aos objectos com que o 
transformador contacta, como sejam os apoios 
em que assenta e os cabos de ligação. Esta 
parte do calor perdido é mínima e pode ser sis- 
tematicamente desprezada nos cálculos. A con- 
dução é também responsável pela transmissão de 
praticamente a totalidade do calor das perdas 
através das paredes da cuba e do permutador 
desde o óleo até à superfície exterior, de onde 
é perdido por convecção e radiação. 

Seja dado um paralelipípedo de pequena espes- 
sura E, cujas bases têm a área A e estão a tem- 
peraturas uniformes cuja diferença é 0), Se o 
paralelipípedo for constituído por um material 
homogéneo de condutividade térmica k, a potên- 
cia calorífica transmitida através dele por condu- 
ção é dada por 

Q=kA (9) (33) 
Nos casos das chapas de aço que constituem as 
cubas e permutadores a condutividade é alta e a 
área A é grande. Sendo limitada a potência Q, o 
quociente ((/E) é pequeno. Como a própria es- 
pessura E é pequena, a diferença de tempera- 
tura O) também o é, podendo em geral ser des- 
prezada nos cálculos (ver a fig. 2). 
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Enquanto a condução se refere em especial à 
transmissão de calor em meios sólidos, a con- 
vecção diz respeito à transmissão entre a super- 
fície dum sólido e um fluido, que neste caso está 
normalmente em movimento. 

Quando um fluido se move em presença de 
um corpo sólido forma-se junto à superfície uma 
camada do fluido em que o movimento é con- 
dicionado localmente pela presença do sólido. 

Dá-se-lhe o nome de camada-limite, para a 
distinguir do fluido não perturbado. Na parte da 
camada-limite que contacta com a superfície o 
movimento é do bem conhecido tipo laminar, 
em que as velocidades locais são constantes no 
tempo e de direcção paralela à superfície. Na 
parte externa o movimento é em geral turbu- 
lento, caracterizado por ter velocidades locais 
variáveis no tempo em módulo e direcção, sendo 
eventualmente oblíquas à superfície. Se esta está 
a uma temperatura diferente da do fluido, troca 
calor com ele. Esse calor transmite-se por con- 
dução através da camada laminar até à zona tur- 
bulenta, onde há um transporte intenso de fluido 
quente que leva consigo o calor que contém. É 
este o mecanismo da convecção, em que a trans- 
missão do calor é consequência dum transporte 
de fluido. 

O movimento pode ser provocado pela im- 
pulsão hidrostática resultante da diferença de 
densidades entre o fluido não perturbado e a ca- 
mada-limite a temperatura diferente, ou por 
qualquer causa puramente mecânica. Conforme 
o caso assim a convecção se denomina natural 
ou forçada. 


Para terminar diremos que é costume expri- 
mir a potência calorífica Q, transmitida por con- 
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